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ABSTRACT 

Optimal regulation of hyperbolic systems in the presence of unknown disturbances is considered. 
Necessary conditions for determining the optimal control that tracks a desired trajectory in the presence 
of the worst possible perturbations are developed. The results also characterize the worst possible 
disturbance that the system will be able to tolerate before any degradation of the system performance. 
Numerical results on the control of a vibrating beam are presented. 

I. INTRODUCTION 

The H 0 o control problem for regulation of dynamical systems in the presence of perturbations has 
been a subject of considerable research in recent years [4,5,6,7,10,16], Although the original formulation 
of the Hoc method was in terms of frequency domain terms, extensions in state space terms leading 
to feedback control using Riccati type equations have been developed [4,7,10,16], For the infinite 
dimensional systems, the II ^ control has started to gain momentum. For a summary of recent results, 
see the survey paper [3]. Like its finite dimensional counterpart, the frequency domain approach [9,15] as 
well as the state space analysis [11,12,] in the presence of both bounded and unbounded perturbations 
has been considered in the literature. The problems pertinent to the II ^ control design are: a) 
input-output stability, b) disturbance decoupling, and c) disturbance attenuation. 

This paper is concerned with disturbance attenuation of hyperbolic systems in the presence of worst 
possible disturbances. We utilize the concepts of optimal control theory [1,8] for infinite dimensional 
systems for deriving the control law for optimum regulation of the system in the presence of worst 
possible disturbances. The method presented in this paper is a generalization of an //^-type method 
developed in [13,14] for finite dimensional systems. The ratio of disturbance energy to the energy of 
the controlled system is used as a measure of performance for disturbance attenuation. We present 
conditions for estimation of the largest perturbation that can be attenuated and the corresponding 
controller to attenuate this perturbation. 

The paper is organized as follows: Section II introduces the control problem. Necessary 
conditions for optimum disturbance attenuation are presented in section III followed by numerical 
results on control of a vibrating beam in section IV. Some concluding remarks are given in section V. 

II. NOTATIONS AND PROBLEM STATEMENT 

We shall use the following notations for abstract function spaces throughout the paper. Let II be a 
Hilbert space, and V a linear subspace of H carrying the structure of a reflexive Banach space with the 
injection VC.H continuous. We identify // with its dual so that VC+HC+V', where V is the topological 
dual of V. 
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Suppose A be a bounded linear self adjoint operator A E C(V, V) satisfying the conditions 

\(A<p,rp)\ < c\\<p\\ v llV’llv, c > 0 , e V 
(A<f,<f)v',v + P\<p\ h > a>0, (3 e R, <p € V 


Consider the hyperbolic system 

d 2 y 

dt 2 


+ Ay — Bu + Cv, t E / = (0 ,T), 

d y, 


2/(0) = 2/o, ^(0) = y' Q 


(2) 


where the operator A is as defined above. The control applied to the system is denoted u£W = L 2 (I , II), 
and B is a bounded linear operator B e C(H). Suppose the system is perturbed by a disturbance 
v E L 2 (I, H) through the operator C E C(H). The initial conditions j/o € V and y' 0 e H are also 
considered to be initial disturbances to the system. 


With this introduction we now pose the control problem: 


Given the perturbed system (2), find the control u E 1*2(1, H) that keeps the state trajectory as 
close as possible to a desired trajectory in the presence of maximum possible additive disturbance 
v G 1*2(1 , H) and maximum possible initial disturbances yo £ V and y' 0 E H. 


For a mathematical formulation of this control problem, we introduce a cost function: 

2*1 In \yo\ 2 dx + \s 2 dx + 3 J Ix nr 2 \v\ 2 dxdt 


J(u,v,yo,y' 0 ) = 


(3) 


291 JixCl \V~ V d \ 2 dx dt + 592 / /x Q |j It ~ yf\ 2 dxdt + \ J /x q n|ti| 2 dxdt 
where si,s 2 ,qi,q 2 ,r\ and r 2 are scalar weighting factors, and y d and yf are desired trajectories respec- 
tively. Then the disturbance rejecting control problem is equivalent to the minimax problem of finding 
a control u and a scalar A* so that 


A ‘ = i|}f sup J(u,v,yo,y'o) (4) 

ueu 

yo^o 

y'o?° 

subject to the dynamics (2). The quantity A* can be interpreted as the disturbance rejection capacity of 
the system. A larger A* implies a better controller in the sense that the system will be able to tolerate 
larger amount perturbations before degradation of the system performance. A small A* means that the 
system is too sensitive to disturbances; despite the effects of the control the state trajectory is not close 
to the desired trajectory even in the presence of a small amount of perturbations. 


We shall assume that a solution of this minimax problem exists. In what follows, we shall derive a 
set of necessary conditions that must be satisfied by the optimal controller. 


III. MAIN RESULTS 


We first give a brief outline of derivation of the main results. The minimax problem introduced 
above is solved in two steps, with the first step being finding the supremum of J over u assuming 
that the perturbations v and yo, y' 0 are known, and the second finding the infimum of J over nonzero 
perturbations. The first step determines the optimal control that regulates the system, and the second 
step characterizes the worst possible perturbation that the controller will be able to attenuate before a 
serious degradation of the system performance. 
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Clearly, the problem of finding the supremum of J over u is equivalent to minimizing the denomi- 
nator of J given in (3) for fixed v, yo, and y' 0 subject to the dynamics (2). This is a well known problem 
in infinite dimensional control theory for hyperbolic systems (see [1,2,8] for details). Theorems 1 and 2 
presented below pertain to this problem. We omit the proofs for brevity. 

Theorem 1. For a given v £ L 2 (I, H), yo £ V, and y' 0 £ H, the system (1) has a unique solution 
y £ L 2 (I,V) n C(I,V), y t £ L 2 (I, H) n C{I,H). Furthermore, the mapping (yo,y’ 0 , u) — ■> (y,y t ) is 
continuous from V x H x L 2 (I,H) — > L 2 (I,V) x L 2 (I,H). ■ 


For convenience of presentation, we introduce two new variables <p\ = y and y> 2 = y t , and rewrite 
the system (2) as a first order equation: 


d<s 

— — I- Ap> = Bu + C v 

dt * 

^(O) = (yo, y'o) = p 


(5) 


where 


0 -/ 

A 0 ’ 



and 


C = 


0 

c ■ 


Using the above notations, we also rewrite the denominator of the cost function (3) as 

Mn)=\j (<p~ <p d ,Q(<p~ V d ))HxH dt + ^ {u, R\u)f{ dt (6) 

where Q = diag (q\,q 2 ) and R\ = n. Let q\ and q 2 be nonnegative, and r\ strictly positive. 

The necessary and sufficient condition that uo £ L 2 {I,H) be optimal in the sense of minimization 
of the cost (6) is that 


J[(uo\ u — uo) > 0 for all u £U 


(7) 


where J[(uo;u - uo) is the Gateaux derivative of J\ at uo £ U in the direction u - uo- This is given in 
the next theorem: 


Theorem 2. Consider the system (5) for fixed additive disturbance v £ L 2 (I,H), fixed initial disturbance 
p £ V x H, and the desired trajectory <p d £ L 2 (I, H ) x L 2 (I, H). Then the optimal control uo £ L 2 (I, H) 
that minimizes the cost (6) is characterized by the solution of the two-point-boundary-value problem: 

^ +A<p + BRf 1 B*i>=Cv, <p(0 ) = p (8) 

-^T + A~ip = Q(<p - <p d ), ffT) = 0 (9) 

The optimal control wo is then given by 

Uo = —Rf l B*\j>. ■ (10) 


At this point we return to the disturbance rejecting control problem introduced earlier. Clearly J(uo) 
is a function of v and p which are yet to be determined. We substitute Ji(wo) into the denominator of 
the cost function (3) leading to 


J, v \(p,Sp)h + \ jfpy^zfyidt 

5 //<¥> - F d , Q(F ~ F d ))H*H dt + \ JfBRpB-f, 0 ) dt 


( 11 ) 
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( 12 ) 


Then the disturbance rejecting control problem is equivalent to finding a scalar A* so that 

A* = inf J(v,p) 

Pf o 

subject to the dynamics (8) - (9). Note that since A* is the optimal solution, we have 0 < A* < A = J(v,p). 
Hence clearly, the following function 

Jl{v,p) = \{p,Sp) H + \ J{ V, R 2 v)h dt - y j{p-p d ,Q{p - <p d ))Hy.H dt - y J (BRp B* rf dt (13) 

is convex and nonnegative, and has a minimum at J 2 = 0. Thus the problem of finding the infimum 
indicated in (12) is equivalent to minimizing (13) subject to the system dynamics (8) - (9). 

By virtue of Theorem 1, it is clear that for any v £ L 2 (I, H ), p £ Vx If, and p d e L 2 (I, II) x L 2 (I, H), 
the equations (8) - (9) have a unique solution p £ L 2 (I,V ) x L 2 {I,H ) and if £ L 2 (I,V ) x L 2 (I,H). In 
addition, the solution has a unique Gateaux derivative satisfying the following theorem: 

Theorem 3. The solution (p, ip) of the two-point-boundary-value problem (8) - (9) corresponding to 
v £ L 2 (I,H) and p £ V x H has a unique Gateaux derivative at every vo £ L 2 (I,II) and po £ V x II 
satisfying 

^+Ap + BRpB*4’ = C(v-v 0 ), p(0 ) = p-po (14) 

-Jf+W-Q$ = 0, 4fT) = 0 (15) 

with p £ L 2 (I, V) x L 2 (I, H) and % £ L 2 (I, V) x L 2 (I, II) ■ 

Necessary conditions of optimality for minimization of (13) are now derived with 
above results and the fact that the Gateaux derivative 

V - VO ,p- po) > 0 

for all v G L 2 (I,H ) and po e V x H, where is the Gateaux derivative at vo,po 
v - vo, p — po • We present the result in the following theorem: 


Theorem 4. The worst additive disturbance vq and the worst initial disturbance po that can be attenuated 
by the optimal control uo are characterized by simultaneous solution of the following equations: 



^ + Ap + BRpB'xp = CR^C’f, 

(17) 


-?jy + = Q(p - </), 

(18) 


~m +A ^~ Qn = 

(19) 


+ At] + BRpB*£ = A *BRpB*if> 

(J t 

(20) 

with the boundary conditions 



^ ~ P ° 4{T) = 0 

£(0) = Spo 

f(T) = 0 

m = o 

(21) 


the help of the 
( 16 ) 

in the direction 
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The worst disturbances are given by 


vo = 

Po = 5-^(0) 

The optimal control that regulates the system in the presence of the worst disturbance is given by 

uq = 

Proof: Taking the Gateaux derivative of the J 2 at r>o G L 2 (I,H) and po £ V x H, we have 

0< J' 2 {voipo',v-vo,p-po) = (Spo,p-po)H + j{R 2 Vo,v-vo)dt-\* j {tp, Q(<p~<f d )) dt- A* Jfy, BRJ l Brp) dt 

The result follows from Theorem (3) and adjoint system (19) - (20). ■ 

It is worthwhile to mention here that equations (17) - (21) represent a two-point-boundary-value 
problem with A* as a parameter which is unknown. The smallest value of A for which this TPBVP, i.e., 
(17) - (21) has a solution is the optimal A* or the disturbance rejection capacity of the system. The 
corresponding control u 0 is then obtained using (23) and the worst disturbance vo and po that can be 
attenuated is given by (22). 

IV. EXAMPLE 

We consider the cantilever beam equation (normalized) 

0 + 0 = d( x )u(t) + h(x)v(t), x e II = (0, 1), t > 0 (24) 

subject to boundary conditions 

2 /( 0 , 0 = 0 , ||( 0,0 = 0 , 0 ( 1,0 = 0 , 0 ( 1,0 = 0 ( 25 ) 

Define the operator A in H = Lziil) by 

A * = 0. D(A) = {* ■.* € - 0, g(0) = 0, 0(1) = 0, 0(1) = 0} 

where H 4 (Q) is the Sobolev space of order four on {}. For V we take V = {4> E // 2 (fi),</>(0) = 

0,f!(0) = 0}. 

We assume that the desired state of the controlled system is the zero state, and that there is no 
initial disturbance. We compute the disturbance rejection capacity of the system using Theorem 4. 
Table I shows that a tighter regulation (i.e., higher Q) is possible only if less disturbance is allowed 
to be attenuated. It is intuitively correct to say that a better regulation of the state trajectory can 
be achieved if the disturbance amplitude is small. Similarly a cheaper control allows more disturbance 
accommodation by the controller as shown in Table II. Stated in a different way, this means that 
attenuation of larger amplitude disturbances will require more control energy. 

TABLE I TABLE II 


( 22 ) 

(23) 


Q 

r\ 

T2 

5 

A* 

1 

1 

1 

10 

0.5522 

10 

1 

1 

10 

0.2130 

20 

1 

1 

10 

0.1521 

50 

1 

1 

10 

0.0922 

100 

1 

1 

10 

0.0618 


Q 

T\ 

T2 

5 

A* 

20 

1 

1 

10 

0.1521 

20 

2 

1 

10 

0.1065 

20 

5 

1 

10 

0.0637 

20 

10 

1 

10 

0.0420 
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V. CONCLUSION 


We present an H^- like control method for hyperbolic systems. Necessary conditions in the form of 
a two-point-boundary-value problem for determining the optimum controller and the worst exogenous 
input that can be attenuated by the optimum controller have been derived. The results are related to 
the //oo control problem in the sense that the H^ norm is given by the inverse of square root of A* 

[14] . The disturbance rejection capacity has been computed for a cantilever beam. Further research 
needs to be done to develop state feedback and output feedback controllers, and to extend the method 
to the infinite horizon control problems. 
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